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By considering f(R) gravity models, the cosmic evolution is modified with respect to the standard 
ACDM scenario. In particular, the thermal history of particles results modified. In this paper, we 
derive the evolution of relics particles (WIMPs) assuming a reliable f(R) cosmological solution and 
taking into account observational constraints. The connection to the PAMELA experiment is also 
discussed. Results are consistent with constraints coming from BICEP2 and PLANCK experiments. 

PACS numbers: 04.50.-h, 98.80.-k, 98.80.Es 


1. INTRODUCTION 

General Relativity (GR) is a self-consistent theory of gravity where space and time are considered as dynamical 
variables and new concepts as black holes and cosmic expansion are introduced. From the cosmological point of view, 
the prediction of cosmic microwave background radiation (CMBR) and the formation of primordial light elements 
maybe represent the greatest success of this theory. However, despite all these fundamental results, GR has not been 
fully investigated at the ultraviolet scales, where strong deviations from the standard Hilbert-Einstein picture emerge, 
and, on the other hand, new ingredients, such as dark matter (DM) and dark energy (DE), are required in order to 
fit the gravitational dynamics at infrared scales [lj. Moreover, self-consistent and comprehensive approaches to deal 
with gravitational interactions at fundamental level (quantum gravity) are still missing. 

These arguments lead to the conclusion that a unitary theory encompassing the gravitational phenomenology at all 
scales is still lacking. In the last years, several alternative or modified theories of gravity have been proposed, also with 
the aim to address the shortcomings related to the Cosmological Standard Model, based on GR. For example, higher 
order curvature invariants than the simple Ricci scalar R allow to get inflationary behaviors, removing the primordial 
singularity, as well as to explain the flatness and horizon problems Q. This approach and, of course, all those related 
to it, are fundamentally motivated by the fact that, at high curvature regimes, further curvature invariants have to be 
considered for constructing self-consistent effective actions in curved spacetime m. In some sense, the introduction 
of higher order terms, depending on the invariants of curvature, are required at high curvature regimes. Of course, 
this is not the final step for building up a quantum gravity theory, but it allows for an effective description that works 
well at least at one-loop level 

All the above motivations strongly suggest that towards ultraviolet regimes (i.e. in high density regimes), GR has 
to be modified by adding further curvature corrections. In the framework of models that extend GR, f(R) gravity 
is certainly one of the favorite candidate since it provides, in a natural way, an almost unified description of DE and 
DM, without invoking exotic sources as DM Q. Moreover, it allows for the unification of the early-time (inflation) 
and the later-time acceleration of the Universe 0,@|. The gravitational action for f(R) gravity is given by 

S =-^ J d 4x V^g f(R) + Smlg^^m], (1-1) 

where S m is the action of the standard matter and « 2 = 87 rG = 87 r/Mp Z , with the Planck mass Mpi ~ 10 19 GeV 
(details and applications are discussed in [8l4llj). This theory can be viewed as a particular case of scalar-tensor 
gravity by dealing with further degrees of freedom as a scalar field (see for details). 
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One of the consequences of dealing with alternative cosmologies, including hence f(R) cosmology, is that the thermal 
history of particles results modified. In fact, one finds that the expansion rates H of the Universe, obtained in modified 
cosmologies, can be written in terms of the expansion rate Hgr obtained in GR, H(T) = A(T)Hgr(T), where the 
factor A(T) encodes the information about the particular model of gravity extending or modifying GR. Usually, the 
factor A(T) is defined in order that the successful predictions of the Big Bang Nucleosynthesis (BBN) are preserved, 
that is A(T) ^ 1 at early time, and A(T) —> 1 before BBN begins (one refers to the pre-BBN epoch since it is not 
directly constrained by cosmological observations). 

The aim of this paper is to explore the possibility of explaining the PAMELA observations, i.e. the measured 
excess in positron flux above ~ lOGeV [h§], in terms of DM thermal relic abundance, by means of f(R) cosmology. 
In fact, since in the framework of the conventional cosmology and particle physics a widely accepted explanation 
for such observations is still lacking, one, in principle, cannot rule out dark matter, as weakly interactive massive 
particles (WIMPs), axions or heavy neutrinos, as a possible solution. It is worth noting that DM interpretation of 
PAMELA data has indeed renewed the interest in alternative cosmologies because, as already noted, these models 
lead in a natural way to a modification of the expansion rate of the Universe jl3| : an enhanced pre-BBN expansion 
can reconcile observed DM cosmic relic abundance with indirect DM detection experiments (such as PAMELA and 
the more recent AMS-02). 

According to the above considerations, we assume that the evolution of the Universe is governed by an f(R) model 
of the form 


f(R) = R + aR n 


( 1 . 2 ) 


This model can be generated, for example, in the framework of Supergravity [14 |. and in a perturbative regime it turns 
out to be a correction of the R + R 2 /M 2 model (extended Starobinsky’s model). Results in Refs. [1J, LL5[ show that 
the Starobinsky model is in good agreement with the BICEP2 data. Moreover, sizable primordial tensor modes can 
be generated in (marginally deformed) models of the form (11.21) . provided 1 < n < 2 jl6| (more specifically in (Tf| it 
is shown that if inflation is driven by f(R) gravity, then a natural form for this function is (I1.2P . where the value of n 
could be related to the microscopic theory dictating the trace-log quantum corrections). The main point is to rewrite 
the f(R) action in the Einstein frame, which implies the appearance of a scalar field. The derivation relies essentially 
on two steps: i ) the introduction of the conformal mode ip = —df/dR and of the real scalar field <j>, related to ip as 
2ip — 1 = £</>, ii) the generation of the kinetic term for <p through the conformal transformation g—> (1 + f,(p 2 )g^- 
Inflation occurs for large values of the scalar field, i.e. <j> £ -1 / 2 . This implies ip 1 or, equivalently, df/dR 1, 

therefore a regime where the i?"-term is dominant. Moreover, the Starobinsky-like inflation model may also emerge 


from dilaton dynamics in brane cosmology scenarios based on string theory 17]. Models based on (11.21) have been 


also studied in the context of bouncing cosmology (see, for example, 18]). It has to be also mentioned that the recent 
analysis by the PLANCK Collaboration [l9j led to the conclusion that I? 2 -inflation ( R 2 /M 2 R) is fully consistent 
with observations @, [20|, . 

At this point is worth a comment about the chameleon mechanism and f(R) gravity [22j. Such a mechanism 
asserts that the Compton wavelength A (typically assumed constant), associated to the characteristic scales obtained 
by adding (pertubative) higher-order terms to the Hilbert-Einstein action, is smaller/larger in those regions where the 
matter density is higher/lower. As a consequence, the theory can be seen as a local effective theory which is valid for 
a certain range of parameters. For f(R) gravity, the chameleon mechanism works because, in the high-energy density 
regions, it reproduces the Newtonian gravitational forces, making the model compatible with the Solar System tests 
gam (see also mi). 

The paper is organized as follows. In Section II we derive the f(R) gravity field equations and solve them in the 
radiation dominated era (well before the BBN onsets), hence from the period of GUT scales to the transition time 
t* >tBBN, when the Universe starts to evolve according to the standard cosmological model. Section III is devoted 
to the study of thermal relics abundance and DM particles required to explain the PAMELA experiment. Conclusions 
are drawn in Section IV. 


2. FIELD EQUATIONS IN f(R) GRAVITY 


The field equations for f(R) gravity follow by varying the action (11.11) with respect to the tensor metric 


= k 2 T™ , G c = f'R, V M V„/' + g^Uf 


df 


where f = and T™ is the energy-momentum tensor for matter. The equation of the trace is 

3 □/' + f'R -2 f = n 2 T m , T m =p~3p. 


( 2 . 1 ) 


( 2 . 2 ) 
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The tensor G£„ satisfies the Bianchi identities, i.e. = 0, so that, for consistency, 

divergenceless too: 

one gets that T™ is 

= 0 . 

(2.3) 

In a (spatially flat) Friedman-Robertson-Walker (FRW) metric 


ds 2 = dt 2 — a 2 (t)[dx 2 + dy 2 + dz 2 ], 

(2.4) 

the nonvanishing components of G c fJ 1 ' are 


G c 0 ° = fR° 0 -\f + ZHf ', 

(2.5) 

G C 1 = f'Ri-tsl + (j' + 4Hf^6i, 

(2.6) 


where we have used □/' = /' + 3 Hf, H = a/a , and the dot stands for d/dt. 

As specified in the Introduction, we work in the regime where the i? n -term dominates ( aR n > R in Eq. (11.21) , with 
1 < n < 2 according to Ref. [l6j]) during the Universe evolution from GUT scales to the transition time t*. The latter 
characterizes the instant in which the Universe passes from the cosmic evolution described by f{R) cosmology to the 
cosmic evolution described by the standard cosmological model (see below). For simplicity, we look for solutions of 
the form a{t) = a^t^. The 0 — 0 field equation and the trace equation give (in the early Universe, hence in the limit 
t -> 0) 


a£l 0<n R n = n 2 p, (2.7) 

aT 0 , n R n = K 2 T m , (2.8) 


where p is the energy density, that in the radiation dominated era reads p 
relativistic degrees of freedom and T is the temperature), meanwhile 


7r2 ff* T 4 
30 


(< 7 * counts the number of 


f^/3 ,n 
T /3,n 


R 


n{f3 + 2n — 3) 


- 1 


= n — 2 — 


2/3-1 
n(n — l)(2n — 1) 


£(2/3-1) 


6/3(1 - 2/3) 
t 2 


3 n(n — 1) 
2/3—1 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


The functions {T^ n , fl 0 n } vs /3 are plotted in Figs. [T]and[51 From (12.71) and (12.81) . it follows that their ratio is given 

, T m Tp ra 

by — = 7^—- 

P *1(3,n 

In what follows we shall consider two cases, T m = 0 and T m ^ 0 (in the first case one gets a relation between /3 
and n, in the second case £ and n can be taken independent): 


T m = 0, i.e. T 0 , n = 0 - In this case, one gets two solutions: 

2n 2 — 3n + 1 


n 


£ i =2> £2 = 

For these solutions, the function assumes the form 


n — 2 


5 n 2 — 8n + 2 n 

^ n= 4(n-l) ’ P= 2 


( 2 . 12 ) 


(2.13) 


and tip 2j „ = 0. The function fl/ 3 1)n is positive for n ^ 1.289. 
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• T m 0 - This possibility may occur, for example, in scenarios where bulk viscosity effects are considered 1 . 
These effects are generated owing to the rapid expansion/compression of fluids, ceasing to be in thermodynamical 
equilibrium. This occurs, in particular, in an expanding Universe, when fluids are out of equilibrium. Typically, 
these processes are so rapid that the system undergoes thermal equilibrium very quickly. However, in the case 





FIG. 2: ra, n vs j3 for different values of n = 1.1,1.5, 1.9. 


i 


It is interesting to note that —— = (1 — 3w) y 0, with p = wp, can be obtained in the case in which, for example, the interactions 
among massless particles are taken into account [31 1. These give rise to a trace anomaly T m oc 3(g) F /,iy F lltJ ^ 0, so that the adiabatic 


index turns out to be w = 


1 -? 


where 


3 


_5_g 4 (N C + jN f ) (fN c - §1V» 

18-rr 2 ( 47 r ) 2 „ , 7 N c Nf 

" + 2 IV 2 - 1 


+ 0 ( 9 5 )- 


At high energies and for typical gauge groups and matter content in the Universe, one finds that the order of magnitude of this quantity 
is given by c ~ 10 —2 — 10 —1 , playing therefore a not negligible role in the early phases of the Universe evolution. Another possibility to 
a have a non vanishing T m is to consider quantum fluctuations of primordial fields [33 . For a FRW Universe one gets T rn = —3fcs/4t 4 , 
where £3 = 1 44 Q 7r ^ (No +317Vi + llA’ 1 / 2 / 2 ), where Ni is the number of quantum fields (for SU(5) model one finds k% ~ 10 —2 ). However, 
this case cannot be used in the present paper, since field equations are not fulfilled. 
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FIG. 3: 7 vs P for different values of n = {1.5, 1.8, 2}. 


in which one considers particle decays of one or more species (see for example [33j]), then a finite time is required 
for driving the system at the equilibrium. The energy-momentum tensor in presence of bulk viscosity term is 


given by [3J, |35| (see also 


T;r = (p+P + IL) uX - (p + n) s ;, 


(2.14) 


with trace T m = —311. Here = (1, 0) is the four-velocity of the fluid ( u 2 = u a ii a = 1), and n the bulk viscous 
pressure. Following [ 35, [H, Q, we concern here the case in which the bulk viscous pressure is proportional to 

p, i.e. n = — 7 p/ 3 . Hence, one gets 7 = ^’ n > 0. In Fig. [3] is plotted 7 vs /? for n = {1.5,1.8, 2}. 

^ */3,n 

Let us now determine the relation between the cosmic time t and the temperature T. From Eq. m one obtains 


where 2 


t = £ 


—y 

Mpi) 


M 


-1 
pi > 


£^[6|/3(l-2/3)|] 1 / 2 


2(1 -n) 


M 


pi 


(2.15) 


(2.16) 


The transition time (temperature) f* (T*) is determined by equating the equation of the evolution in f(R) cosmology, 
Eq. (l2T7ll . with that one in GR, i.e. 

a£lp, n R n {U:) = Hq R (U) . 


One gets 

t* = [4&fy,n[6|0(2/9 - 1 )|]"]^tj M ~} ? 
The expression of the transition temperature T* is then given by 


T* = 


15 

167r 3 p* 


—[ 4a ^’ ra ] 4( lDi— M P i 
[6(|/3(1- 2/3)|)] 


(2.17) 


(2.18) 


2 We have introduced the absolute value on the quantity (3(2(3 — 1) because of reality of t* for all (3. Although the field equations remain 
the same under the change R —► —R, the quantity (3(2/3 — 1) change sign in order that the R n is always well defined for all f3. 
















that allows to recast the relation (12.171) in the form 
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Moreover, notice that 

UTl _ I 15 
Mpi V 167T 3 g* 

At the end we recover that the expansion rate of the Universe in f(R) cosmology can be written as 

H(T) = A(T)H GR (T), A(T) = 2V3p(^-) , u=--2 

\1* J n 

where the factor A(T) is the so called enhancement factor. 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


3. RELIC ABUNDANCE AND WIMP PARTICLES 


In this Section we study the thermal relic abundance in f(R) gravity. As mentioned in the Introduction, alternative 
cosmologies indeed predict modified thermal histories of relic particles that occur during the pre BBN epoch. When 
the expansion rate of the Universe changes, as compared to that one derived in the framework of GR, thermal 
relics decouple with larger relic abundance, a scenario that might have its imprint on relic WIMPs (weakly interacting 
massive particles). The interest about these particles as DM follows from the fact that WIMPs in chemical equilibrium 
have the same abundance of cold DM. These studies are also motivated by astrophysical results obtained for cosmic 
ray electron and positrons 0 . l40l - l42l |. antiprotons [43j| . and y-rays OEH- As mentioned, particular attention is 
devoted to PAMELA experiment that observed a excess of positron fraction at energies greater than lOGeV 0 ]. 
Besides the astrophysical interpretation of this phenomenon proposed in Refs. 0 ], there is also the possibility that 
the raise of the positron fraction could be ascribed to DM annihilation into leptons [47], EH . In this last case, a large 
value of (tJannv) is required. More specifically, PAMELA and ATIC data require a cross section of the order or larger 
than (<Jannv) ~ 10 -26 cm 3 sec -1 . Such a value is also necessary in order that thermal relics have the observed DM 
density (see also mi)- 

The characteristics of the Universe expansion, such as the composition and/or the expansion rate, affect the relic 
density of WIMPs (and more generally of other DM candidates) as well as their velocity distributions before structure 
formation. According to the standard cosmology and particle physics, the calculation of the relic density of particles 
relies on the assumption that the period of the Universe dominated by radiation began before the main production of 
relics and that the entropy of matter is conserved during this epoch and the successive one. However, any contribution 
to the energy density (in matter and geometrical sector) modifies the Hubble expansion rate, and, as a consequence, 
the value of the relic density. Investigations along these lines have been performed in different cosmological scenarios 
01111 . 


The general analysis that accounts for the enhancement of the expansion rates in alternative cosmology has been 
performed in Ref. 0 (see also mi). The expansion rate H is written in the form H = A(T)H G r , where the 
function 3 A(T) = rj(T/Tf) v encodes, through the free parameters {v,r]}, a particular cosmological model, Tf is the 
temperature at which the WIMPs DM freezes-out in the standard cosmology, Tf ~ lOGeV (Tf, in general, varies by 
varying the DM mass m x ). The parameter v labels cosmological models: v = 2 in Randall-Sundrum type II brane 


cosmology [53| , v = 1 in the kination models 


v = 0 in cosmologies with an overall boost of the Hubble expansion 


rate [13|, v = —1 in scalar-tensor cosmology [55f. In our f(R) model, we have v = 2 jn — 2, so that — 1 ^ v ^ 0 for 


1 < n < 2. 


3 In 0, the enhancement function A(T) is parameterized as 


A{T) = 


1 + r ?(^)" tanh2 7W £ 


for T > T bbn 
for T < T BB jy 


(3.1) 


where T BB n ~ IMeV. This form of A(T) allows to avoid conflicts with BBN. In the regime T 3 > T BB ,y. the function 13.11 behaviors 
as A(T) ~ rj . 
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TABLE I: In this Table are reported some estimations of a for fixed values of the transition temperature T* = (1 — 10 2 )MeV. 
The expressions of /3 and are given in (12.131) . 


n 

T» (MeV) 

a 

1.3 

1 

10 14 GeV -0 ' 6 


10 2 

10 12 GeV -0 ' 6 

2 

1 

10 44 GeV -2 


10 2 

10 36 GeV -2 


For our estimations, we shall refer to the analysis performed in jl3| . where the conditions for which modified 
cosmologies can explain both the PAMELA data and the abundance of relics particles, without violating the constraints 
provided by astrophysical observations, have been studied. More specifically the analysis concerns the DM annihilation 
cross section (a ann v) vs DM mass in the interval [lOGeV-lOTeV] (these annihilation channels are typical of several DM 
particles, such as lightest SUSY or Kaluza-Klein particles), for different annihilation channels (e + e - , t + t - , 

W + W~, bb) and for different DM density profiles (Via Lactea and Aquarius DM distributions). The study is performed 
by numerically solving the Boltzmann equation for the number density of thermal relic 4 Y, taking into account for 
the modifications related to the expansion rate H = A(T)Hgr■ Fixing v = {—1,0,1,2} and {a a n n v) ~ 2.1 x 10 -26 
cm 3 sec -1 , one determines the values of the parameter 77 vs m x = [lOGeV, lOTeV], required to infer the correct relic 
abundance of DM particles fl x h 2 = D x /i 2 |^^ P = 0.1131 ± 0.0034 [56). 

The analysis in [{ 3 } shows that the values of the parameter 77 necessary to explain the PAMELA data (in particular 
for the case of DM annihilation into e + e - ), together with tt x h 2 ~ 0.11, are (see Figg. 11, 12, and 15 of [13| )) 

77 ^ 0.1 for 777 x > 10 2 GeV. (3.2) 

More precisely, for DM masses in the range [10 2 —10 4 ]GeV, the allowed region for the parameter 77 is 0.1 ^ 77 < 10 3 —10 6 , 
where the upper bounds on 77 vary for the different cosmological models labelled by v. 


Applications to f(R) cosmology 

According to the above results, we rewrite the factor A(T) (see Eq. (12.211) 1 in the following form 


A(T) = 

V = 

v = 



(3.3) 


The transition temperature T* is fixed for values greater than the free-out temperature Tf. Therefore we set T* = 
(1 4- 10 2 )Tbbn- From (12.181) we get 


a = 


15 

167T 3 g* 


"- 1 [6|/3(l-2/3)|]-» /MpA 4( " -1) . 

4 n 0 ,n \ r» / 


4 


The relic abundance is given by D x ^ 2 = m x^o y o ^ w j iere p c — SH^Mp^Sn is the critical density of the Universe, so is the present value 
of the entropy density, and Yo is the present value of the WIMP abundance for comoving volume f55ll 


1 

Lb 


1 




g X {x){CannV) ^ 
Vg*(x)A(x)x 2 


X = 


m x 

T 


Here Yf is the value of the WIMP abundance for comoving volume at the freeze-out, {g x (T), g*(T)} counts the effective number of 

Mpim x (a ann v)f 


degrees of freedom at temperature T, and Xf = In 


0.0038<7x" 


, which is computed for non relativistic DM particles. 


' A(xf)y/xfg*(xf)_ 

Notice that in f(R) cosmology one has x = qxH , where q = n/2/3. This factor does not alter the Boltzmann equation since q = 1 for 
(3 = n/2 and it is of the order q ~ G( 1 — 2), for the values (3 ~ 1 and 1.3 < n < 2 here used. We can therefore safety use results of Ref. 

0 . 
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FIG. 4: r/ vs n for /? = n/2 and transition temperatures T = {1,10,10 2 }MeV. T/ = lOGeV is the freeze-out temperature, while 
r; = 0.1 is the lower bound on r/, see Eq. (IQIl . 



The order of magnitudes of a are reported in Table I for /? and fgiven in (12.1311 . 

The function rj vs n is plotted in Fig. [I] for £lp >n and /3 = n/2 given in (12.131) . corresponding to T m = 0, and in 
Fig. [5] for flf 3 ;n given in (12.91) . corresponding to T m ^ 0. In both cases, the parameter ?y assumes values of the order 
0(0.1 — 1), so that the mass of WIMPs particles is of the order 10 2 GeV. Notice finally that the enhancement factor 
(1331) increases for larger values of /3, hence for a super-accelerated expansion of the early Universe. 


4. CONCLUSIONS 

In this paper we have studied the problem of thermal relic particles in f(R) cosmology. We focus on a power 
law model of the form f(R) = R + aR n , which is consistent with recent PLANCK Collaboration and BICEP2 data 
constraints: n = 2 corresponds to Starobinsky’s model, while the so called marginally deformed model (n ^ 1) 
produces sizable primordial tensor modes provided the exponent n falls down in the range n G [1,2]. As we have 
showed, if the cosmic evolution of the early Universe is described by modified field equations, as provided indeed 
by f(R) gravity, then the expansion rate gets modified by a factor A(T) ( H(T ) = A{T)Hgr(T ), see Eq. (12.211) 1. 
This quantity essentially weights how much the expansion rate of the Universe in f(R) cosmology deviates from the 
expansion rate derived in the standard cosmology, and affects, in turn, the production of relic particles (thermal relics 
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decouple with larger relic abundances). As a consequence, the latter is obtained for larger annihilation cross section, 
and therefore also the indirect detection rates get enhanced. This effect may have its imprint on supersymmetric 
candidates for DM. 

For a power law scale factor, solutions of the modified field equations, and parameterizing the enhancement factor 
f T\ v 

as A(T) = rj[ — , we find that the f(R) model is consistent with PAMELA data (for DM annihilation into lepton 

\ T fJ 

channel e + e _ ), and the abundance of relic DM f l x h 2 ~ 0.11, provided that 77 ~ 0(0.1 — 1) (and — 1 ^ v < —0.46). 
According to (13.21) . the corresponding WIMPs masses are m x > 10 2 GeV. 

Finally, it is worth noticing that the analysis here performed relies on the model in which the form of f(R) is a 
power-law expansion of the scalar curvature R, as well as on power law solution of scale factor. In general, it would 
be interesting to consider other curvature invariants, such as the Riemann, the Ricci and the Gauss-Bonnet tensors, 
and their derivatives, which could play a relevant dynamical role for the evolution of relic particles, as well as to seek 
a more general solution of the field equations for the f(R) model here discussed. 
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